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ABSTRACT 

We calculate the synchrotron self-Compton emission from internal shocks occurring in relativistic winds as 
a source of gamma-ray bursts, with allowance for self-absorption. For plausible model parameters most pulses 
within a Gamma-Ray Burst (GRB) are optically thick to synchrotron self-absorption at the frequency at which 
most electrons radiate. Up-scattering of photon number spectra harder than j/° (such as the self-absorbed emission) 
yields inverse Compton photon number spectra that are flat, therefore our model has the potential of explaining 
the low-energy indices harder than v^^^^ (the optically thin synchrotron limit) that have been observed in some 
bursts. The optical counterparts of the model bursts are sufficiently bright to be detected by such experiments as 
LOTIS, unless the magnetic field is well below equipartition. 

Subject headings: gamma-rays: bursts - methods: numerical - radiation mechanisms: non-thermal 



1. INTRODUCTION 

In the framework of internal shocks occurring in highly rel- 
ativistic winds, the up-scattering of synchrotron photons was 
taken into consideration in the calculation GRB spectra by Pap- 
athanassiou & Meszaros (1996) and Pilla & Loeb (1998). Both 
groups found a reasonable qualitative agreement between mod- 
eled and observed spectra. As shown in §p^, for electron in- 
jection fractions of order unity and reasonable dissipation ef- 
ficiencies, synchrotron emission is liable to peak below the 
spectral peak observed in GRBs at hundreds of keVs (Band et 
al. 1993), even for equipartition magnetic fields. For this rea- 
son we consider a GRB model where the synchrotron peak is 
in/around the optical range, and the 7-photons arise from in- 
verse Compton (IC) scatterings (Meszaros & Rees 1994). In 
this case synchrotron self-absorption (SSA) can be important 
at the synchrotron peak, yielding up-scattered spectra that are 
harder below 100 keV than an optically thin synchrotron emis- 
sion peaking at hard X-raysFI. 

In this work we obtain numerically the synchrotron self- 
Compton emission from internal shocks through simulations 
of the wind dynamical evolution and the emission released in 
each collision. We also investigate the brightness of the optical 
flashes radiated by internal shocks (Meszaros & Rees 1997). 
IC scattering of synchrotron seed photons has been previously 
considered as a mechanism of emission of high energy pho- 
tons in GRB afterglows (e.g. Meszaros & Rees 1994, Sari, 
Narayan, & Piran 1996, Waxman 1997, Chiang & Dermer 
1999, Panaitescu & Kumar 2000) and in blazars (e.g. Boettcher, 
Mause, & Schlickeiser 1997, Mastichiadis & Kirk 1997, Urry 
1998, Chiaberge & Ghiselhni 1999). 

2. DESCRIPTION OF THE MODEL 

2.1. Wind Dynamics 

The relativistic wind is approximated as a sequence of dis- 
crete shells of uniform density, ejected by the GRB source with 

^ If a photospheric component is present, then thermal and pair-Comptonized 



various initial Lorentz factors and masses. We assume that the 
time it takes the GRB engine to become unstable and eject a 
shell is proportional to the energy of that shell. This implies 
a constant wind luminosity L on average throughout the entire 
wind ejection duration t^,. We consider that the shell Lorentz 
factors r have a log-normal distribution (Beloborodov 2000), 
r — 1 = (Fq — 1) exp(^a;) with x Gaussian distributed, and A 
of order unity. 

The radii at which the collisions between pairs of neighbor- 
ing shells take place are calculated from the ejection kinemat- 
ics, taking into account the progressive merging, the adiabatic 
losses between collisions, and the radiative losses. The Lorentz 
factor and internal energy of the shocked fluid are calculated 
from energy and momentum conservation in each shell colli- 
sion. The shock jump conditions (Blandford & McKee 1976) 
determine the speeds of each shock, which give: i) the shocked 
fluid energy density u' (used for calculating the magnetic field 
B and minimum electron Lorentz factor 7^, see §2.2); ii) the 



shell shock-crossing time t'^ (which is the duration of the in- 
jection of elec tron s, determining the cooling electron Lorentz 
factor 7c, see §2.3); and Hi) the post-shock shell thickness A' 
(used for calculating the adiabatic cooling timescale and the 
radiative efficiency). Between collisions the co-moving frame 
thickness increases at the sound speed. 

2.2. Emission of Radiation 

The synchrotron self-Compton emission from each shocked 
shell is calculated from the Lorentz factor, magnetic field, num- 
ber of electrons (given by the shell mass) and the electron dis- 
tribution in each shocked shell. The turbulent magnetic field 
strength B and the typical electron Lorentz factor ji resulting 
from shock-acceleration are parameterized through their frac- 
tional energy ^ and e relative to the internal energy: 
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emission could also produce hard spectra (Meszaros & Rees 2000) 
1 



where u' and nj. are the internal energy and particle densityQ 
respectively, and p is the index of the power-law distribution 
of injected electrons: Mi{'^) oc 7^^ for 7 > 7^. Denoting 
by 6 the laboratory frame internal-to-kinetic energy ratio in 
the shocked fluid (i.e. the dissipation efficiency), one can write 
u' ~ 6n'm„c^, hence 



is i;(7) - (7mec2)/(Piy + Pl^, where P^ - YP'^y is the 
electron IC radiating power and P' = (l/67r)crThC7^i3^ is 
the synchrotron power, thus 
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with An = 10^"yl and using p = 3 for simplicity. 

The density tiJ. can be calculated from the shell mass m, 
thickness A', and radius r, using m = Aiir^mpn'^A' . For an- 
alytical calculations one can approximate the typical mass m 
of a shell radiating a pulse which arrives at observer time T 
as a fraction T/T[, of the total wind mass M, where Ti, is the 
burst duration, which is approximately equal to the duration 
tw of the wind ejection. One obtains m = LT/Tc^, where 
L = E/t^, is the wind luminosity, E and F being the wind to- 
tal energy average Lorentz factor of the wind. The collisional 
radius r can be derived from the spread in the photon arrival 
time 5 = r/T'^c due to the spherical shape of the emitting sur- 
face, which is roughly the pulse duration^. The shell thick- 
ness can be approximatedn by A' = r/T. Therefore we obtain 



= LT/{'iTTmpe'T^5^). Then equation ^ yields 
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If the emission at i^i is not self-absorbed, then the peak fre- 
quency Vsy of the synchrotron power vFJ)''^ is at the character- 
istic synchrotron frequency i^i = (ei?7j?F)/(4mec) of the 7^- 

(ic) 

electrons, and the peak of the IC power i/Fi) is at u, 



From equations (E) and (51), one obtains that 
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t'^ being the time elapsed since the beginning of injection of rel- 
ativistic electrons, i.e. the shell shock-crossing time, which is 
calculated numerically from the shock dynamics. 

From equation (H) the cooling electron Lorentz factor jc de- 
fined by 4(7,) = t^ is 7e - Ci/(y + 1). For 7, < 7, 
the 7i-electrons cool faster than their injection timescale. In 
this fast electron cooling regime the resulting electron distri- 
bution has a low-energy tail J\f{'j) oc 7^^ for ^c < 1 < I'h 
while 7V(7) oc 7~(p+i) above 7^. In the slow electron cool- 
ing regime, characterized by 7^ < 7^, the electron distribu- 
tion is the injected one, M{j) oc j^^, for 7i < 7 < 7c, and 
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above 7c, due to cooling. As physical pa- 



rameters vary among shocked shells, the GRB pulses can be in 
different electron cooUng regimes. 

2.4. Compton Parameter and Self-Absorption Break 

We set the electron distribution as described in §2.3. The 
Compton parameter is given by 
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where t^ is the optical depth to electron scattering and Ta{'^) 
is the optical thickness to SSA at the synchrotron characteris- 
tic frequency for 7-electrons, and 7^ is defined by 77i(7a) = 1. 
The first integral in the right-hand side of equation (Q) takes into 
account that only those photons emitted within a region (around 
a given electron) of optical depth unity are up-scattered before 
being absorbed. Table 1 lists the values of the Y parameter for 
all possible orderings of 7a, 7c, and 7^, assuming that these 
breaks are sufficiently apart from each other, and ignoring mul- 
tiplicative factors of order unity. 

For a power-law distribution of electrons it can be shown 
that the optical thickness Ta to SSA at the characteristic syn- 
chrotron frequency for electrons at the top of the distribution, 
i.e. at 7p = min(7c, 7^), is 
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For slow electron cooling (7p = 7^) one obtains 
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Equations (0) and (ph show that the peak frequencies of 
the optically thin synchrotron and of the IC emissions de- 
pend strongly on parameters which may vary substantially from 
pulse to pulse. Equation (0) also shows that for plausible model 
parameters the synchrotron peak hvsy cannot be as high as 100 
keV for pulses longer than about (5 = 0.1 s. 

Due to the collisions among shells, the amplitude of the 
Lorentz factor fluctuations in the wind is reduced, diminishing 
the dissipation efficiency 6. Equation (ph shows that the peak 
of the IC emission depends strongly on and thus should have 
a general trend of decreasing in time. This could be the rea- 
son for the overall spectral softening seen in GRBs (e.g. Ford 
et al. 1995). 

2.3. Electron Distribution 

The electron distribution N{'^) at the end of shell energiza- 
tion is determined by the Afi{'y) initially injected at shock, by 
the electron radiative cooling (through synchrotron and IC ra- 
diation) and by the absorption of the synchrotron photons. The 
radiative cooling timescale of electrons with Lorentz factor 7 

^ Primed quantities are in tlie co-moving frame. The values used in the numerical calculations are those when the shock has swept up the entire shell. 
^ The true pulse duration is larger, due to contributions from the shell shock-crossing and the cooling time (radiative and adiabatic) of the electrons radiating at the 
observing frequency. These are taken into account numerically and influence the light-curves and instantaneous spectra, but not the burst-integrated spectra. 
* This is a good approximation for a freely expanding shell; shock-compression and shell merging yield smaller shell thicknesses. 



If electrons are cooling fast ("/p = 7c), then Tp is even larger. 
Equation (^ shows that is quite likely that the synchrotron 
emission from the typical, 7p-electron is self-absorbed, thus the 
synchrotron spectrum Fi, peaks at the synchrotron frequency Va 
corresponding to ja- 

The shell optical thickness to SSA at any frequency 1/ (or 
corresponding 7) can be calculated from Tp given in equa- 
tion (ph by using Ta{h') oc v^^/"^ oc ^^^^/"^ for 7 < 7p and 



Ta{v) oc v '•'?+^-'/^ oc 7 '•'?+^) for 7p < 7, where g is the 
index of the electron distribution around 7, i.e. q = 2 for 
7c < 7 < 7i, 9 = P for 7i < 7 < 7c, and q = p + 1 for 
max(7i,7c < 7. 

To take into account that for jc ^ 7a the synchrotron cool- 
ing is basically suppressed by self-absorption, we calculate jc 
from 7c = Ci/[Y + x(7c/7a)], where x is a function satisfy- 
ing lim^j^o x{x) = 0, x(l) = 1 - e"\ and lim^j^oo xi^) = 1- 
Note that the Y parameter and the self-absorption 7^ depend 
on the electron distribution, which is at its turn determined by 
the electron cooling through synchrotron (possibly reduced by 
self-absorption) and IC. Thus the equations for 7^, Y (see Table 
1), and 7c are coupled. For given parameters (Ci, C2,Te,7i,p) 
these equations are solved numerically for each of the cases 
listed in Table 1 and only the self-consistent solution (i.e. the 
one that satisfies the assumed ordering of 7a, 7c, and •ji) is re- 
tained for the calculation of the synchrotron and IC emissions. 

2.5. Synchrotron and Inverse Compton Spectra 

The synchrotron spectrum fI)"^' is approximated as a se- 
quence of four power-laws with breaks at I'a, Vi, and Vc- The 
slope of the spectrum between break depends on the ordering 



of these breaks: Fij 



(sy) 



Any) 



oc v 



oc v^''^ for niin(i/c7 i^j) < v < Va, 
^ for u < min ( z^a, t^c, i^i)' F^ oc v^^^ for 



v„ < V < niin(i'c, t'i), -Fl;^^ oc v ^'^ for max (i^q, i^c) < 
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The up-scattered spectrum has breaks and slopes that are de- 
termined by those of the synchrotron spectrum F^' and of 
the electron distribution Af{'j), thus the IC spectrum is more 
complex than that of the synchrotron emission. Figure I shows 
the logarithmic derivatives (i.e. slopes) of the IC emission for 
all possible orderings of the synchrotron breaks, obtained by 
integrating numerically over F^^' and N'{'f) the up-scattered 
emission per electron given in equation (2.48) of Blumenthal 
& Gould (1970). Note that the IC spectrum resulting from 

?isy} 



Fy oc v^, by a typical electron with 7 ~ min(7c, 7^) cannot 

iic) 

be harder than Fl) oc v. 

Numerically, we calculate the IC spectrum by approximating 
it as a sequence of power-laws, with breaks at the frequencies 
identified in Figure 1, and we ignore higher order IC scatter- 
ings. For the model parameters we shall consider, scatterings 
of third order or higher are suppressed as they occur in ex- 
treme Klein-Nishina regime. However this is not always true 
for second IC scatterings. Therefore our calculations ignore a 
very high energy (above 1 GeV) component, and may overesti- 
mate the brightness of the synchrotron and first IC components 
(but not their ratio) as some electron energy would be "drained" 
through a second up-scattering. 

3. SIMULATED BURSTS 

The burst spectrum is calculated by summing the synchrotron 
and IC spectra from all collisions. Figure 2 shows the time- 
integrated spectrum $1, = J Fi,{T)dT for a wind consisting of 
13 shells, for plausible model parameters (L, t^,; e, <^), and for 



a log-normal distribution of F. The burst has two major struc- 
tures whose spectra are shown separately. The first one (0 — 5 s) 
consists of two pulses whose IC emissions peak at 10 keV and 
2 MeV and has a rather soft spectrum in the BATSE range. The 
second structure (15 — 20 s) is a single pulse, self-absorbed at 
the synchrotron peak, with the IC emission peaking at 200 ke V 
and having a hard spectrum below 40 keV. The other six colli- 
sions occurring in the wind considered here are very dim in the 
10 keV-1 MeV range either because their radiative efficiency 
is low or because their synchrotron and IC emissions peak far 
from this photon energy range. 

Other sets of parameters may also lead to hard low en- 
ergy slopes resulting from up-scattering of self-absorbed syn- 
chrotron emission. However, since the physics of the model 
does not confine the peak frequencies IC emission of the bright 
pulses within a decade, the spectra resulting from the addition 
of many pulses do not exhibit in general a steep slope below 
and around 10 keV. Therefore our model for the GRB emission 
can explain the hard low energy spectra observed by Preece et 
al. (1998) only for bursts with a modest number of pulses. 

Figure 3 shows the dependence of the i?-band magnitude of 
the burst optical counterpart on the most important model pa- 
rameters. The wind parameters used in these calculations pro- 
duce GRBs of average intensity and with spectra compatible 
with the observations. Notice that most of the optical counter- 
parts are above the reported sensitivity of the up-graded LOTIS 
experiment (Williams et al. 2000). The exception is the case of 
a wind with a very low magnetic field parameter £_. 

The non-detection of such optical counterparts can be used 
to infer lower limits on the IC parameter, which sets the relative 
intensity of the high (7-ray) and low (optical) frequency emis- 
sion. It is straightforward to show that a "typical" GRB lasting 
10 s, with a vFiy peak frequency around 100 keV and fluence of 
5 X 10^^ ergcm^'^, has an i?-band magnitude 



i? = 8.8 + 2.5 (1 - /?) log X + 2.5 log Fa 
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where X is the ratio of the peak frequency of the synchrotron 



up-scatterings of photons below Va, where F^"^' oc i/^/^ or power vF^, 
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to that of observations (vf; 
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^ 4.7 X 1014 Hz), 

and /3 is the slope of F^""' , i.e. F^ oc v^^. Equation ( pX]| ) 
shows that, if the synchrotron peak frequency is in or not too 
far from the optical domain, as is the case for the winds whose 
optical light-curves are shown in Figure 3, then the average R- 
magnitude of the optical flash is dimmer than Rlotis ^ 13 
(WilHams et al. 1999) if F > 5000. 

For fast electron cooling and in the Va < Vi case, i.e. for 
t he fi rst two cases^ given in Table 1, the electron distribution of 
§2.3 yields Y ~ jcJiTe which, in the ^ <C e limit, leads to 
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where u'g = B^ 



(11) 
TT is the magnetic energy density, and A' ~ 
ct'^ was used. The above condition for LOTIS non-detection 
(Y > 5000) impHes ^ < 10"* £-.5. Therefore the optical flash 
of a "typical" burst would be dimmer than about _R = 13 if the 
magnetic field is several orders of magnitude below equiparti- 
tior 



^ These cases are encountered in a fraction of the simulated pulses for the parameters given in Figure 3, nevertheless the increase of Y with decreasing g is a 
general feature exhibited by the simulated bursts. 

® We note that if only a small fraction of the injected electrons were to acquire the fractional energy e and if the magnetic field is close to equipartition (^ ^ 1), 
such that the synchrotron emission would peak around 100 keV, the optical flashes would also be dimmer than the LOTIS limit. 



4. CONCLUSIONS 



The analytical treatment presented in §2.2 and the numeri- 
cal spectrum displayed in Figure 2 show that IC up-scattering 
of synchrotron emission from internal shocks in relativistic un- 
stable winds can produce GRBs with break energies and low- 
and high-energy indices that are typical of real GRBs (Preece 
et al. 2000). The peak frequency of the IC emission is strongly 
dependent on some of the model parameters (see eq. [|]| for the 
first two cases listed in Table 1). For a range of plausible pa- 
rameters this peak falls within the BATSE observing window. 

The "harder than synchrotron" low-energy spectra that have 
been reported in a significant fraction of bursts (Preece et 
al. 1998) can be explained by the Compton up-scattering of syn- 
chrotron spectra that are self-absorbed. Even though the syn- 
chrotron spectrum can be as hard as Fu oc v^'^ below the ab- 
sorption break, the up-scattered spectrum cannot be harder than 
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therefore the self-absorbed self-Compton model 



has a limiting value aic < (the low energy index of the photon 
spectrum). This is larger by da = 2/3 than the limiting value 
for an optically thin synchrotron spectrum, leaving only a cou- 
ple of bursts in Figure 2 of Preece et al. (1998) with as exceed- 
ing aic by more than Icr. We emphasize that self-absorption of 



the synchrotron emission is not a guarantee that the resulting 
model spectrum is as a hard as a,c — 0, as the addition of the 
IC emission of many pulses with various IC peak frequencies 
may result in a flatter burst spectrum, and that, in general, hard 
low energy spectra can be obtained only for bursts with a few 
pulses. 

The prompt optical to gamma-ray emission ratio from inter- 
nal shocks depends on the fractional energy in the magnetic 
field - a poorly known parameter. This is independent and ad- 
ditional to a possible optical flash from a reverse component 
of a subsequent external shock (Meszaros & Rees 1997, Sari 
& Piran 1999). A possible explanation for the non-detection 
of optical counterparts down to i? ~ 13 by LOTIS is that the 
magnetic field in internal shocks is several orders of magnitude 
weaker than the equipartition value, corresponding to an aver- 
age Compton parameter around or larger than 1000. 
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TABLE 1. 
The Compton Parameter (Y), self- absorption Lorentz factor (7a) and peak of 

THE IC l^Fi, spectrum (l^ic) FOR VARIOUS POSSIBLE ORDERINGS OF 7a, 7c, AND -fi 
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Fig. L — Slopes of IC spectra resulting from up-scattering of a broken power-law (four segments) synchrotron spectrum by a broken power-law distribution of 
electrons, for various possible combinations of synchrotron breaks: Ua (absorption), Ui (injection), and i/c (cooling). Horizontal lines indicate the slopes of the 
synchrotron spectrum, which may be asymptotic values for the IC spectrum slope. Note that the IC emission below the lowest energy break indicated for each 

(ic) 

spectrum (red arrows), which results from up-scattering of self-absorbed synchrotron emission, is .F^; oc i/. Blue arrows/characters indicate the IC break frequency 

(ic) 

closest to the peak of the i/F^) emission. Black an'ows indicate less conspicuous breaks or IC frequencies that separate two asymptotic values of the IC spectral 
slope. A steep injected electron distribution with p = 4 was used so that the IC breaks are easier to identify. IC spectra were calculated in the co-moving frame and 
blue-shifted by a factor 100. 
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Fig. 2. — Time-integrated spectrum for a smmlated burst produced by a wind consisting initially of 13 shells, whose ejection Lorentz factors has a log-normal 
distrihutipn with To = 320 and A = 1 (see ^M. Other wind parameters are: t„ = 10 s, L = 3 X 10^^ ergs^^, £ = 0.15, § = 10"^, p = 3, redshift z = I 
(see §E3. The upper panel shows separately the spectra of the two structures exhibited by the light-curve, identified by the time ranges when they are seen, as well 
as the overall burst spectrum. Note the steepness of the 15 — 20 s spectrum below 40 keV peak (lower panel), due to up-scattering of self-absorbed synchrotron 
emission. The average dissipation efficiency per collision is 6 = 13%. Approximately 29% of the wind energy is eventually dissipated, out of which ~ 4% is stored 
in electrons. About 50% of the electron energy is lost adiabatically, therefore the bolometric radiative efficiency of the burst is only 2.0%. The efficiency at which the 
wind converts its kinetic energy into 20 keV-1 MeV photons is 1.2%, due to that only 60% of the total emission falls in this energy range. 
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Fig. 3. — i?-band magnitude of burst optical counterparts for few inodel parameters, bined in 2 s time intervals. The reference set is tu, = 10 s, L = 10^^ ergs~^, 
e = 0.15, 5 = 10~^, p = 3, redshift 2 = 1 and 100 ejected shells. Fq = 320 and A = 1 for F log-normal distributed, while for the uniform deviate F is randomly 
distributed between 50 and 10"^. The bursts last for about (1 + z)tw = 20 s, exhibit 20-35 peaks, have spectral features (low and high energy spectral slope, break 
frequency) consistent to those found by Preece et al. 2000 in real bursts, 20keV-l MeV fluences in the range (2 — 8) X 10^^ crgcm^^, and efficiencies of radiating 
the wind total energy in the BATSE range between 0.3% and 1%. 



